Magnetoresistance of mesoscopic granular ferromagnets 
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We have performed magnetoresistance (MR) measurements of granular ferromagnets having lat- 
eral dimensions smaller than 0.5 fim and containing a small number of grains (down to about 100). 
Compared to macroscopic samples, these granular samples exhibit unusually large saturation fields 
and MR amplitudes. In addition, the evolution of the magnetoresistance curve as the intergrain 
distance decreases is qualitatively different than that of large samples. We discuss these results and 
suggest that they reflect a transition from percolation to quasi single-channel dominated transport. 

PACS numbers: 73.23.-b; 73.40.Rw; 75.50.Cc 



I. INTRODUCTION 



Insulating granular ferromagnets, i.e. systems of mag- 
netic grains imbedded in a non-magnetic insulating ma- 
trix, exhibit negative magnetoresistance (MR) curves. 
At zero magnetic field (H=0) the resistance is maximal. 
For applied ± H the resistance decreases until, for large 
enough field, it reaches saturation. This behavior was ob- 
served in a variety of granular samples 0, 0> El El IE IE 
and has been ascribed to spin dependent tunneling be- 
tween randomly oriented magnetic moments of the grains 
[Elil- Applying a magnetic field aligns these moments, 
causing the tunneling resistance to decrease. The magne- 
toresistance is thus determined by the relative magnetic 
orientation of pairs of grains. The MR amplitude (de- 
fined as AR/R = R ^ H ~*^Q^ R ^ ) °f a pair of grains i and 
j is given by HE]: 
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where P is the electron polarization and a is the angle 
between the moment orientations at zero magnetic field. 
In the presence of magnetic field, H, the magnetic ori- 
entation of each grain is governed by two factors. The 
first is the magnetization easy axis of the grain due to its 
anisotropy, and the second is the external magnetic field. 
The total magnetic energy of a grain per unit volume is 
given by 
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where fj,Q is magnetic permeability in vacuum, v is the 
anisotropy coefficient (of the order of 1), M s is the satu- 
ration magnetization, is the angle between the magnetic 
moment orientation and the easy axis of magnetization 
and <p is the angle between the magnetic moment and 
the applied magnetic field. Eq. 2 is composed of two 
energies. The left component, Wi, represents the energy 
due to orientation out of the easy axis and the right com- 
ponent, Wh is due to alignment of the moment with an 



external field. The magnetic orientation of each grain 
(and therefore the MR amplitude) at zero temperature 
is determined by minimizing the energy of equation 2. It 
should be noted that ultrasmall grains at finite temper- 
atures undergo a superparamagnetic transition, thus the 
thermal energy, kgT, may dominate and overcome the 
effect of the magnetic energy. 

When dealing with transport properties through a 
granular insulator, one has to take into account that not 
all the grains participate in the conductance processes 
since it is a strongly disordered system. It has long been 
realized that a percolation treatment is the proper way 
to deal with such a system as it provides much insight 
into the physics of the conductivity [l2l Hii , fl4| . In this 
approach each pair of grains i and j is represented by a 
resistor with resistance Rij inversely proportional to the 
hopping probability between the grains and given by |l5| : 
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where is the distance between the grains, £ is the lo- 
calization length representing the decay of the electronic 
wavefunction in the insulator and Eij is the energy dif- 
ference between the electronic states. For metallic grains 
at temperatures of the order of a few K, E.^ is related 
to the charging energies of the grains which depend on 
the grain diameters. In general, smaller grains give rise 
to larger EijS. The granular system can be mapped by 
a resistor network containing series and parallel current 
paths. Because the network contains a wide distribution 
of resistances (due to the exponential factors in Eq. 3) 
the transport is percolative and governed by a set of criti- 
cal resistors. These act as "red bonds" of the percolation 
network and their conductivity determines the transport 
properties of the entire system. Hence, the scale of in- 
homogeniety is the percolation radius, Lc, which can be 
viewed as the average distance between critical resistors. 

As the lateral dimension of the sample is reduced be- 
low Lc, the nature of the transport is expected to change 
dramatically. A percolation network is no longer a suit- 
able way to treat the system. Instead, a single current 
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path, or even a single critical resistor, is expected to dom- 
inate the transport, giving rise to mesoscopic effects and 
large sample to sample fluctuations. 

In this paper we study the MR curves of granular fer- 
romagnets with sizes smaller than Lc (~ 0.5/xm) and 
compare them to the properties of macroscopic samples. 
We find that reducing the size causes a dramatic increase 
in both MR amplitudes and saturation fields. In addi- 
tion, the evolution of the MR curve as a function of aver- 
age inter-grain separation is different than that of large 
samples. We discuss these results and attribute them to 
a crossover from percolation transport to a single dom- 
inating current trajectory, as the sample size is reduced 
below the percolation radius. 

I. EXPERIMENTAL 

The samples described in this paper were discon- 
tinuous Ni films prepared by "quench condensation" 
0, El EH- i- e - evaporation of thin, discontinuous 
films on substrates that are kept at cold temperatures 
and under UHV conditions as described elsewhere 0, . 
This techniques enables one to perform electric and mag- 
netic measurements during the sample growth. Thus, one 
can study the magneto-transport properties of a single 
granular Ni sample as a function of film thickness. We 
note that the thickness barely changes during the experi- 
ment. Adding ~ A to a film of nominal thickness 30 A is 
sufficient to reduce the sample resistance by a few orders 
of magnitude Q ■ This demonstrates that the decrease in 
resistance occurs due to increasing of inter-grain coupling 
while the grain sizes remain practically constant during 
the sample growth process. Thus the quench condensa- 
tion method enables one to study the magneto-transport 
of granular ferromagnets as a function of the mean inter- 
grain distance without thermally cycling the sample or 
exposing it to atmospheric conditions. 

In order to study small sized granular samples we 
combined photo- lithography and atomic force microscope 
(AFM) fabrication methods with quench-condensation. 
First we prepared Ni wires having width of a few /im on 
a Si02 substrate by conventional lithography. Next we 
used an AFM tip (with the z feedback loop disabled) to 
cut the wire in two, thus creating two close electrodes. 
The nano-space formed between the electrodes defined 
the geometry of the measured sample. The properties 
of the nano-space were determined by the quality of the 
AFM's tip, the force that was applied on the wire and 
the smoothness of the Ni. A typical electrode configu- 
ration is shown in the inset of fig. 1. Once a desired 
geometry was achieved, the substrate was placed on a 
quench-condensation probe and a granular film was evap- 
orated into the gap. Using this technique we were able 
to prepare samples with sizes as small as a few tens of 
nm. Since the grain diameters are of the order of 10-15 
nm [20I |2~0 ] these samples contain about 100 grains. We 
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FIG. 1: AR/R for a 3mm * 3mm sample (top panel) and 
a lOOnm * lOOnm sample (bottom panel). The insert is an 
AFM image of the mesoscopic electrode template. A gap with 
lateral dimensions of lOOnm is cut from a Ni wire having width 
of 3(im. 

name all samples with lateral dimensions smaller than 
0.5/im "mesoscopic samples", while larger samples are 
named "macroscopic samples" 

I. RESULTS 

Magnetoresistance (MR) curves of two samples, one 
macroscopic and one mesoscopic, having sheet resistance 
of 2 MQ, (nominal thickness of ~ 21 A), are depicted in 
Fig. 1. The magnetic field in these experiments was ap- 
plied perpendicular to the film plane. The MR amplitude 
of the macroscopic sample is 2% and the saturation field, 
Hs, is ~ 0.5T. These values are typical of all our macro- 
scopic 2D granular samples and are similar to those ob- 
tained by other groups using different preparation meth- 
ods for fabrication of insulating granular ferromagnets 
II II II A saturation field of 0.5T can indeed be 
expected since it is close to H = 4ttM s , M s being the 
saturation magnetization, which is the field required to 
align a thin Ni film perpendicular to the substarte against 
the shape anisotropy. 

Unlike macroscopic samples, the mesoscopic samples 
exhibited large sample to sample variations of both MR 
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FIG. 2: AR/R for different samples having sheet resistance 
of 2AIQ as a function of sample lateral size 

amplitudes and saturation fields. Hs is always larger 
than the typical value of 0.5T and varies between IT to 
fields higher than 6T (the largest available field). AR/R 
also fluctuates substantially from sample to sample, how- 
ever, the average MR amplitude sharply increases with 
decreasing sample size as seen in Fig. 2. 

As material is added to the sample and the sheet re- 
sistance, R, decreases .several trends are observed in the 
macroscopic samples @, Q . The saturation field, H s , re- 
mains constant at a value of ~ 0.5T throughout the entire 
sample growth process (Fig. 3c). The MR amplitude, on 
the other hand, decreases monotonically and smoothly, 
until, for R < 0.5kU, it is suppressed altogether (Fig. 
3a). This can be expected since adding material causes 
coalescence of grains. For thick enough layers the sample 
is simply a continuous Ni film in which no spin dependent 
tunneling resistance is expected. 

Another observed trend is the splitting of the MR peak. 
For high resistance the MR curve is centered at H=0 as 
seen in Fig. 1. Upon decreasing resistance, the curve 
splits into two peaks and a hysteresis develops in the MR 
curve H ■ The coercivity, H c grows as a function of de- 
creasing R and reaches 0.25T for the lowest resistance in 
which MR is measurable (Fig. 3b). This was attributed 
to the increase of the effective grain size as the resistance 
decreases, thus giving rise to a transition from a non- 
hysteretic superparamagnetic sample to a ferromagnetic 
sample 6]. 

The mesoscopic samples show qualitatively different 
results. The MR amplitude does not decrease throughout 
most of the sample growth process. On the contrary, 
AR/R appears to increase initially as material is added 
to the sample (Fig. 3a). Only for resistances smaller 
than 50fcf2 a sharp decrease in the magnetoresistance is 
observed. A similar effect is seen for the saturation field, 
Hs- It remains constant for most of the growth and 
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FIG. 3: The evolution of AR/R (top panel), the coerciv- 
ity He (middle panel), and the saturation field Hs (bottom 
panel), as a function of resistance obtained as material is 
added to the sample. Full squares are for a macroscopic sam- 
ple and open symbols are for a number of mesoscopic samples. 



decreases sharply (apparently approaching the bulk value 
of H = AttMs ~ 0.5T) when the resistance drops below 
a few kfl as seen in Fig. 3c. 

Mesoscopic samples differ from macroscopic ones in the 
evolution of He as well. The mesoscopic systems show 
no development of hysteresis as a function of resistance. 
Some of our samples have no hysteresis for the initial 
deposition stages, and no hysteresis is observed even for 
the lowest measured resistance. Other samples exhibit 
a two peak MR curve even for the initial evaporation 
stages. In these samples the initial coercivity does not 
increase as material is added to the system (see Fig. 3b) . 

I. DISCUSSION 

In order to understand these results we note that our 
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FIG. 4: Simulation of the current trajectories in a mesoscopic 
sample containing 180 grains, showing a single domination 
current trajectory. Two leads are connected to the sample 
on the left and right edges, with a small d.c. voltage applied 
between them. The line width is proportional to the current 
flow magnitude. 



mesoscopic systems are smaller than the percolation ra- 
dius, Lc- We have performed computer simulations of 
the MR of our granular Ni containing a small number of 
grains j^. For systems with less than about 500 grains 
we find that the transport is governed by a single dom- 
inating channel. The current trajectories branch out, 
forming a complex percolating network, only for larger 
samples. This corresponds to Lc of about 25 grains 
across, which, in our samples, is equivalent to a sample 
size of about 0.5 /im. For smaller samples the current is 
forced to flow through a single dominating chain of grains 
(see Fig. 4). The mesoscopic nature of these sample ac- 
count for the large sample to sample variations of the 
MR curves in the small samples. We suggest that our 
results can be understood if we assume that the same 
chain of grains dominates the transport throughout the 
sample growth, from a configuration of electrically iso- 
lated grains until the sample is close to the metallic state. 
Adding material does not significantly alter this current 
path. The resistance decreases because the inter-grain 
distance is reduced but the MR amplitude does not re- 
duce until the grains coalesce, the film becomes contin- 
uous and tunnelling MR is no longer expected. This is 
very different to the situation in large samples in which 
adding material opens up new channels and modifies the 
percolation network. 

One consequence of the above model is that the elec- 
tric current in small samples may flow through relatively 
small grains. In macroscopic samples, hopping through 
very small grains is energetically unfavorable due to a 
large contribution to Eij in Eq. 3. The current can 



bypass very small grains and choose larger grains to con- 
struct the transport network. This was shown in refer- 
ence @ where the coercivity extracted from MR mea- 
surements was found to be larger than that extracted 
from magnetization measurements demonstrating that 
the percolation network is constructed of grains larger 
than the average. In small samples the situation is dif- 
ferent since the transport may be forced to take place 
through smaller grains. The fact that only few grains 
are present between the electrodes forces the conduct- 
ing electrons to hop to small grains with large charging 
energy. This may explain the large saturation fields ob- 
served in our mesoscopic samples. The magnetization of 
grains having diameters smaller than lOnm was shown to 
consist of ferromagnetically aligned core spins and a spin- 
glass-like surface layer |23j . Canting of the surface spins 
introduces magnetic stiffness of the grain moment since 
aligning them requires an extremely large external mag- 
netic field. Another reason for large H$ in small grains 
is the affect of thermal fluctuations on the magnetic mo- 
ment orientation. As shown in Eq. 2, the magnetic 
energy is proportional to the grain volume. For small, 
superparamagnetic grains, the energy due to the exter- 
nal H, Wh, has to overcome the thermal energy rather 
than the energy due to the easy axis, Wj. Therefore, 
the smaller the particle the larger the affects of surface 
moments and temperature, giving rise to large H$. As 
the granular sample size is reduced, the average grain 
size participating in the transport decreases. This may 
account for the large saturation fields in our mesoscopic 
samples. Indeed, we observe a clear correlation between 
large saturation fields and small hysteresis in the MR 
curve (see Fig. 5) reflecting the role played by small 
grains in causing large Hs- 

The behavior of the He can also be interpreted apply- 
ing the above considerations. The crossover from non- 
hysteretic to hysteretic MR curves in large samples is 
attributed to an increase in the average grain size in the 
current network. Adding material to a granular array 
causes clustering of the transport network, adding par- 
allel trajectories which allow the current to flow through 
larger grains. He thus increases as the resistance of the 
sample is decreased. In mesoscopic samples this process 
does not occur since the current path is relatively fixed 
throughout the sample growth process. Hence, the grains 
that participate in the transport for high R are expected 
to dominate at low R as well (though the inter-grain tun- 
neling rate is larger) . This accounts for the constant He 
observed in the mesoscopic samples. 

A more subtle issue is that of the MR magnitude. The 
average MR of a random array of Ni grains (according to 
Eq. 1) is expected to be (~ j^pn — 1)- Our simulations 
[22| show that for our Ni films this should produce an av- 
erage AR/R value of about 10% Our macroscopic 
samples exhibit significantly smaller values. In addition, 
Fig. 2 shows that the smaller the sample the larger is 
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FIG. 5: The coercive field, He, determined from the field 
position of the MR peak, as a function of the saturation field, 
Hs, for different small samples (solid squares) and a large 
sample (empty circle). The solid line is a guide to the eye. 



the average MR value. The reason for this is not clear. 
Naively, one could expect that a denser network would 
emphasize the importance of trajectories with small mis- 
match of magnetic moment orientation. Our simulations 
show, however, that the average MR amplitude does not 
depend on the number of parallel channels. This discrep- 
ancy between experiments and theory impels us to sug- 
gest that the MR in the granular ferromagnets is affected 
by magneto-static dipole-dipole interactions which cause 
alignment of the magnetic moments of the grains even in 
the absence of magnetic field. Indeed, we have observed 
signs for magnetic interactions in different granular Ni 
systems |2lj . If ferromagnetic correlations are significant 
in the granular sample, they may lead to a reduction of 
the measured AR/R relative to the theoretical expecta- 
tions. In this case Eq. 2, representing the total magnetic 
energy of a grain, i, should contain an additional factor 
due to magnetic interaction with a neighbor grain, j: 



Wgg 
Vi 



-fxoMg-J-cos(aH) 
% 



(4) 



One can expect that such interactions will reduce the 
randomness of the initial magnetization orientations thus 
suppressing the MR amplitude. We note, however, that 
if the grains are small, the thermal energy can be strong 
enough to smear out the interaction effects in a similar 
manner to the effect on the easy axis as grains become 
superparamagnetic. Therefore, for smaller grains, the 
MR amplitude can be larger than that of larger grains. 
In our mesoscopic samples, in which the average grain 
size is smaller, the interactions will play a smaller role 
and the MR magnitude may be larger, as indeed ob- 
served in the experiments. This does not depend on the 
quench condensed sample resistance since the same set 
of grains dominate the transport throughout the sample 
growth until the film approaches the metallic phase. In 
macroscopic samples, adding material increases the av- 
erage grain size participating in the transport as noted 
above. The importance of the magnetic interactions (rel- 
ative to the thermal energy, KbT) thus increases with 
decreasing resistance. This trend manifests itself in the 
reduction of AR/R with film growth as demonstrated in 
Fig. 3a. 

In summery we point out that reducing the size of 
any disordered system, so that it enters the mesoscopic 
regime, is always accompanied by rich phenomena. The 
case where the sample is ferromagnetic introduces novel 
effects in addition to the usual sample to sample fluctu- 
ations. These include large MR amplitudes, large sat- 
uration fields and unique coercivity behavior. We sug- 
gest that the experimental findings are indicative of a 
transition from percolativc transport to an effective ID 
conductance. In the magnetic case, this enables the de- 
tection of the presence of ultra-small grains which are 
"invisible" in large samples. Assuming magnetic inter- 
grain correlations enables us to account for the different 
MR amplitudes measured in macroscopic and mesoscopic 
samples. Obviously, A more detailed theoretical treat- 
ment of the issues discussed in this paper is needed, this 
is the subject of an ongoing study. 
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dation (grant number 326/02). 



where Vi and Vj are the grain volumes and dij is the 
distance between their centers. This energy should be 
compared to Wi which dominates at H=0. The ratio be- 
tween the two energies for two identical grains in contact, 
having radius a and volume V, is given by: 



, Wgg _ 2 V cos(a H ) 
Wi ' ~ v{2aY sin 2 {9) 



C 



(5) 



where C is a constant of the order of unity. Thus, the 
contribution of the interactions to the magnetic energy 
is of the same order of that arising from the easy axis. 
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